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The paper presents an adaptive control technique for a damaged large transport air- 
craft subject to unknown atmospheric disturbances such as wind gust or turbulence. It is 
assumed that the damage results in vertical tail loss with no rudder authority, which is 
replaced with a differential thrust input. The proposed technique uses the adaptive pre- 
diction based control design in conjunction with the time scale separation principle, based 
on the singular perturbation theory. The application of later is necessitated by the fact 
that the engine response to a throttle command is substantially slow that the angular rate 
dynamics of the aircraft. It is shown that this control technique guarantees the stability 
of the closed-loop system and the tracking of a given reference model. The simulation 
example shows the benefits of the approach. 

I. Introduction 

The survivability of damaged aircraft is one of the main research topics in the aviation safety research. 
The catastrophic loss of control of American Airlines Flight 587 illustrates the degree of hazards we can 
expect due to structural failures of airframe components. 2 Nevertheless, there are many stories in aviation 
history of aircraft surviving after suffering major structural damage to their airframes. The recent example 
of survived aircraft is the DHL A300 cargo flight, which completely lost the hydraulics after it was struck 
by a missile after departing from Baghdad. 1 Pilots were able to stabilize and to land the aircraft using 
differential thrust. 

Using differential thrust as an emergency substitute for failed control surfaces has been a research topic 
after the development of Propulsion-controlled aircraft (PCA) system at NASA Dryden Research Center, 
which was first evaluated on a piloted B-720 simulation. 9 The research on the PCA system continued at 
NASA Dryden and Ames Research Centers in simulations as well as in actual flight test for different flight 
platforms. 4 The objective of the PCA system is to provide a necessary thrust command for each engine for 
emergency flight control in response to pilots flightpath input, assuming that each engine can be controlled 
individually. 

The majority of results in this held are based on known model of the vertical tail loss aircraft, and 
derive thrust control command using conventional control methods (see for example References 3,7,10,21 ). In 
Reference 14 an adaptive control scheme is used to control the aircraft in the presence of unknown actuator 
failure, when differential trust is used to generate necessary moment. In all this cases the engine dynamics 
are ignored, despite the fact that the engine response has a significant lag, especially in law trust levels. 

In this paper we consider a propulsion control problem for a transport aircraft with vertical tail loss, 
taking into account the aircraft’s engine dynamics. It is assumed that the damage alters the stability and 
control derivatives and mass properties of the aircraft. In addition, the aircraft dynamic model involves 
modeling uncertainties that account for aerodynamic forces and moments, and external disturbances. Since 
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the aircraft is missing the rudder, it is assumed that the necessary yaw moment for the directional control 
is generated through the differential thrust of the engines, the dynamics of which is know to be much slower 
than the rotational dynamics of the aircraft. Moreover, it is assumed that the engine dynamics can be 
modeled as an unknown stable system with throttle input and thrust output, which can be measured (or 
very accurate estimates are available). Our adaptive control approach uses several ideas from L\ adaptive 
control, 5 , 6 singular perturbation method 13 and adaptive time scale separation 12 to derive aileron, elevator 
and right and left engine throttle commands to track a given reference command. It is based on the state 
prediction model based control design, which makes it possible to introduce a low-pass filter in the control 
channel to suppress the high frequency oscillations that may be generated by adaptive laws. 5 For the roll 
and pitch channels this approach generates acceptable control signals. For the yaw channel the necessary 
differential throttle input has to be generated as an output feedback signal of a cascaded system, comprised 
of the yaw dynamics and engine dynamics. Therefore, it is tempting to use the the control architecture from 
the Reference 6 that allows to introduce a simple model for the engine dynamics, thus simplifying the output 
feedback adaptive design. However, this approach guarantees good tracking when the model is sufficiently 
fast, which can not be justified for the slow engine dynamics. Therefore, we propose to use the singular 
perturbations method in adaptive settings, which is based on the prediction model as well. 12 The required 
differential thrust is computed to control the prediction model, and the differential throttle command is 
computed via the time scale separation, which is proved to be valid for the prediction models constructed 
for the yaw rate dynamics and engine dynamics. Afterwards, the closeness of the models to corresponding 
dynamics are proved. Also, it is proved that the closed-loop system is stable and the tracking of the reference 
commands can be achieved in steady state as well as in the transient. 

The rest of the paper is organized as follows. In Section II we give the problem characteristics and 
assumptions, and formulate the control objective. The engine dynamics is discussed in Section III. The 
corresponding error dynamics are derived in Section IV. Section V presents the control design steps, and the 
stability analysis is presented in Section ??. Simulation examples are presented in Section VI, and the paper 
is concluded by Section VII. Throughout the paper bold symbols are used for vectors and small letters for 
scalars. 


II. Problem Formulation 

Let the position vectors of the center of gravity Oq of undamaged aircraft be ro in some inertial frame 
Ff. When the damage occurs, the center of gravity shifts to a new position Cft, given by the vector r\ in the 
same inertial frame. It is assumed that the shift is constant and is described by the vector p = r\ — r$. As 
a consequence of this shift a coupling appears between the aircraft’s linear and angular accelerations. The 
corresponding dynamic equations are derived in 16 and have the following form: 

F — m{v + u x v) + mu x (u x p) + moj x p (1) 

M = Jqu: + uj x Jqu> + mu> x (p x v) + mv x (u ; x p) + mp x v — p x W 

These equations are written with respect to a body frame Fb, which is attached to the nominal center 
of gravity Oo, with the longitudinal axis through the aircraft’s nose. In these equations, v is the velocity 
vector of the point Oo and is the aircraft’s angular rate vector, both expressed in the frame Fb, F 
and M represent the sum of external forces and moments in the same frame, and W is the gravity force. 
m = m* + A?n is the unknown mass of the damaged aircraft, where the mass change Am from the known 
nominal mass m* is negative. Jo is the unknown moment of inertia of the damaged aircraft, computed in the 
frame Fb- Substituting v from the first equation in (1) into the second one, after some algebra the angular 
rate equation can be written as 

M — pxF + pxW — (Jo 4- mpp T — ?nl 3 x 3 p T p)A + w x J 0 u + mp T uj(p x us) . (2) 

We notice that J = Jo + mpp T — ml 3 X 3 p T p is the moment of inertia of the damaged aircraft computed in 
the parallel to Fb frame attached to the new center of gravity 0\. We decompose the external force F into 
the sum of the aerodynamic force Fa, gravity force W and the collective thrust T generated by left and 
right engines. 


F = Fa + W + T. 


(3) 



Similarly, we decompose the external moment M into the sum of the aerodynamic moment M a, the moment 
of the gravity force W and the moment generated by the differential thrust AT. The aerodynamic moment 
Ma is further decomposed into the sum of the moment generated by the control surfaces Mg and the 
moment generated by the rest of the aircraft Mg. 

M A = M b + M s + pxW + axAT. (4) 

where a denotes the radius vector of the center of the right engine thrust vector. Substituting the force and 
moment equations into (2), the angular rate dynamics can be expressed in the following form 

Juj = - uj x J q uj + mp T uj(p x uj) — p x (T - W) + Mg + a x AT + M b — p x Fa , (5) 

Next, we substitute uj into first equation in (1) and obtain a force equation in the following form 

F 

v = uj x v — ujuj t p — J -1 (— uj x Jquj + mp T uj{p x uj) 

m 

- p x (T — W) + Mg + a x AT + M b - p x Fa) x p (6) 

The complete dynamics of the damaged aircraft includes also the orientation equations and navigation 
equations. However, in this paper we are not interested in navigation problems as well as the yaw orientation 
equation. That is, the damaged aircraft’s dynamics are given by the equations (5), (6) and 

c t> = p + q sin (j) tan 9 + r cos <f> tan 9 

9 = q cos cj> — r sin (j> (7) 

We assume that the thrust vector is aligned with the longitudinal axis, that is T = [T 0 0] T and 
AT = [AT 0 0] T , where T and AT are the magnitudes of the collective thrust and differential thrust 
respectively. This quantities are subject to the following constraints 

0 < T < T max (8) 

-\t< at < ir , 

where T max is the given maximum collective thrust that can be generated by engines. The second inequality 
in (8) implies that the one sided thrust is allowed during the maneuvers. 

It is assumed that the damage of the aircraft results in the vertical tail loss. This in turn results in the 
directional instability of the aircraft with no rudder authority. The objective is to control the aircraft using 
the available control surfaces and the differential thrust. 

III. Engine Dynamics 

In general, typical aircraft engine dynamics are highly nonlinear and of distributed character. The 
response depends not only on the throttle input and other combustion parameters, but also on the thrust 
level. Figure 1 displays some examples of engine response data to a full throttle input at different altitudes 
with initial thrust levels corresponding to trim conditions for a generic transport aircraft model. 8 

It also shows that for the low thrust level the initial delay in response is much larger than on the high 
thrust level, and the response is faster at law altitudes. The forms of the engine response data suggest that it 
is possible to approximate the engine dynamics by a stable linear system with time delay, which afterwards 
can be transformed into a high order linear time invariant system using Pade approximations. Figure 2 shows 
the forth order linear approximation of the engine response at low trust level, which corresponds to a high 
altitude cruise. It can be seen that the high order linear time invariant approximation describes the engine 
response with sufficient accuracy. Therefore, for the purpose of this paper it is assumed that the thrust 
generated by engines can be described as the output of a slow but stable minimum phase linear system with 
known relative degree and with a single throttle input. 

Since the left and right engines are assumed to be identical, the collective thrust and the differential 
thrust are described by the same dynamic model 

e~ 1 x e (t) = A e x e (t) + b e (u e (t) + d e (t)) 
y e (t) = c e x e (t) . 


(9) 




Figure 1. Engine response data for different altitudes and initial thrust level 



Figure 2. Engine response data and 4th order linear approximation response at altitude 30000 ft 


where x e (t ) is the model’s state, A e , B e , C e are matrices of compatible dimensions with A e being Hurwitz, 
e is a small parameter indicating that the engines dynamics are slow in the real time scale, and d e (t) is a 
bounded external disturbance. For the collective thrust generation y e (t) = T(f), u e (t) = ctr(i), where <5-r(i) 
is the collective throttle input, and for the differential thrust generation y e (t) = A T(t), u e (t) = S& r(t), 
where <5 at(£) is the differential throttle input. The throttle inputs to the left and right engines are defined 
as 


5t l — — [<5t + ^at] (10) 

$t r = 2 — ^ AT ] > 

IV. Error Dynamics 

With the vertical tail loss the aircraft is still symmetrical with respect to the xz plane, that is the center of 
gravity shift has zero y coordinate and J yz = 0. Therefore, it is possible to decouple the longitudinal and the 
lateral-directional dynamics. However, coupling remains in the control inputs via the thrust generation. The 
point will be clarified shortly. Keeping this in mind, we linearize the aircraft’s dynamics about a trim point 
that corresponds to a damaged state, and write the longitudinal and lateral-directional dynamics separately 


x L {t) = A L x L (t ) + B l [A L u L (t) + f L (x L (t),x D (t)) + d L (t )] 


( 11 ) 


where the longitudinal state xl and the control input ul are defined as 


i r i 

xl = AV A a q A 6 , wi = T 6 e , 

the control allocation matrix Al £ R 4x2 is assumed two ne known, matrices matrix Al £ R 4x4 and Al £ 
R 2x2 are unknown, where is diagonal positive definite, and represents the unknown control effectiveness, 
the nominal value of Al is the identity matrix. The continuous and locally Lipschitz function f l (xl,xd) 
represents modeling uncertainties and the bounded time function di(f) is added to account for external 
disturbances such as the turbulence or wind gust. In the above equation, 

r i T 

xl — f3 r p (j> 

is the lateral-directional state that satisfies the equation 

x D (t) = A D x D (t) + B D [A D u D (t) + f D (x L (t),x D (t)) + d D (t)\ , (12) 

where ud = [d 0 AT] t is the lateral-directional control input, the matrices Ad, Bd , A d and functions 
/ d (xl,xd), doit) have the same meaning as in the case of longitudinal dynamics. 

Remark 1 . In the above equations the collective thrust T and the differential thrust AT are considered as 
independent input signals that are generated through the engines, the required thrusts of which are defined 
according to the equations 

Tieft = T + ^AT (13) 

T right =T-^AT, (14) 

where T and AT are considered commands generated from the perspective of controlling the longitudinal and 
lateral- directional dynamics respectively. In this case the collective thrust remains intact while the engines 
undergo changes to generate the required differential thrust. Therefore, the decoupling of the dynamics is 
complete. However, as we will see in the simulations, this way of generating the differential thrust is not 
always the best solution, since for the vertical tail loss aircraft the directional stability is a priority and all 
possible control capabilities must be used for this purposes. The alternative is to control the engines separately 
to achieve maximal differential thrust at each operating point in the expense of altering the desired collective 
thrust, which cannot be considered as an independent control signal any more. The excess collective thrust 
will affect the longitudinal motion, which can be compensated for after the directional stability is archived. 

We construct the reference models for the longitudinal and lateral-directional dynamics that meet the 
desired performance characteristics as follows. Let Al = Af + A Al and Ad = Af, + A Ad, where Af and 
A° n are known matrices such that the pairs (Al, Bff) and (Al, Bff) are controllable, and A Al and A Ad are 
the uncertainties that satisfy the following model matching assumption. 

Assumption 1 . For Hurwitz matrices Alui and Ad m there exist constant matrices Kf x £ R 2x4 , K* Lr £ 
R 2x2 , K* Dx £ R 2x4 and K* Dr £ R 2x2 such that 

B l K* Lx = AA l B d I<* Dx = A Ad 

B l K* Lt = B Lm B D K* Dr = B Dm ■ (15) 

Remark 2 . The true knowledge of the matrices Kf x , Kf r , Kf, x and Kfjr is not required, only their 

existence is assumed. □ 

The reference models are obtained from the corresponding dynamics in the absence of uncertainties and 
disturbances by means of the linear controls 

u ° L (t) = - K Lx x L(t)+K° Lr r L (t) 

u °D(t) = -KD X x D (t)+K° Dr r D (f), (16) 



where the gains iC£ x , and A'£, x . are chosen such the the matrices A^m = A° L — BlK® x and Aom = <~ 
BdKq x are Hurwitz. The reference models for the longitudinal dynamics and lateral directional dynamic 
have form 


A l(^) 

&Dm(t') = AD m XDm(t) A Bj^ m V jj(f) , (17) 

where the r^lt) and ro{t) are external command signals. 

We take the adaptive augmentation approach in the control law design. That is the control signals are 
represented as 

UL(t) = u° L {t) + u a L d {t) 

u D {t ) = u° D (t) , (18) 

where vP L (t) and u° D (t) baseline control signals in (16), and uf d (t), u a p{t) are adaptive augmentations to be 
designed later. Denoting the longitudinal tracking error by e^(f) = x L (t) — XL m (t), the error dynamics can 
be written in the 

e L (t) = A Lm e L (t) A B l [A L (u L (t) - u° L (t )) + K Lx x L (t) A K Lr r L (t ) + f L (x L (t), x D (t)) A d L (t)\ , (19) 

where K Lx = (I - A l)K° Lx + K* Lx and K Lr = (I - h L )K° Lr + K* Lr . Similarly, the lateral-directional tracking 
error dynamics can be written as 

e D (t) = A Dm e D (t ) A B d [A D (u D (t) - u° D (t)) A K Dx x D (t) +K Dr r D (t) A f D (x L (t), x D (t)) A d D (t)\ , (20) 

where e D (t) = x D (t) - x Dm (t), K Dx = (I - A D )K° Dx + K* Dx and K Dr = (I - A D )K° Dr + K* Dr . 

Remark 3 . We notice that the error dynamics (19) and (20) are in the form considered earlier in the 
literature and the suitable adaptive control design techniques guaranteeing tracking both in the steady state 
and in transient can be found in References. 20, 22 Next step would be to design suitable laws for the throttle 
inputs in order to generate the required collective and differential thrust signals designed in the previous step. 
Whereas for the longitudinal design the direct application of any design approach would perform satisfactory, 
for the lateral- directional design it leads to a throttle command that involves high, gain terms. 19 The reason 
is that, the engine dynamics are much faster than the fungoid mode and much slower than the rotational 
dynamics controlled by the rudder deflection. One of the remedies for the situation is the principle of the 
time scale separation, which is based on the Tikhonov’s theorem. 13 This theorem requires some exponential 
stability properties for the cascaded systems under the consideration. However, no exponential stability of the 
error system (20) can be establish without the parameter convergence, which cannot be guaranteed in general, 
especially when the unknown nonlinearities are involved. This prevents the direct application of for the 
time scale separation. As an intermediate step we follow the approach from Reference 12 and design suitable 
prediction models for which the time scale separation theorem can be applied. Since the longitudinal dynamics 
do not require two time scale separation, we drop them from further consideration, assuming that the suitable 
longitudinal design is already available, and the longitudinal variables reside in the suitable bounded regions. 
That is to design control laws for the lateral- directional dynamics all longitudinal variables are assumed to 
be bounded. 

Since f d (xl, xp) is continuous, we approximate it on a compact set f l XL x fl XD by any known approx- 
imation method. Here, we use normalized radial bases functions 11 to approximate f(v,u>) by a linear in 
parameter neural network on the compact set TI Xl x TL Xd . That is we write 

f D (x L , x D ) = W T ip(x L , x D ) + e(x L , x D ) , (21) 

where W £ K JVx2 is the unknown optimal weight matrix, which minimize the norm difference || f d (xl, x d) — 
W T ip(xl,xd) ||, (p(xl,xd) is the vector of basis functions and e(xL,xo) is function reconstruction error that 
is uniformly bounded by a constant e* . We notice that from the above considerations the relationship xl £ 
Vl XL is satisfied, and xjj £ fl XD will be guaranteed by the control design. Substituting the approximation 
(21) into the error dynamics (20) we obtain 

e(t.) = A m e(t) + B [Aw(f) + K x x(t.) + K r r(t) + W T ip{t.) + cr(t)] , 


( 22 ) 



where for the notational simplicity we drop the subscript " D” and introduce notations u(t) = U£,(t) — u° D (t), 
ip(t) = ip(xL(t),XE>(t )) , cr(t) = e(xL(t),XD(t)) + doit), where cr(t) is bounded componentwise by some 
positive constant vector er* as long as the relationship (xL{t),XD(t)) £ fl XL x fl XD holds. Also, since A rn is 
a Hurwitz matrix there exists a symmetric positive definite matrix P that satisfies the Lyapunov equation 

A T m P + PA m = -Q (23) 

for some symmetric positive definite matrix Q. 

V. Control Design 

In this section we construct the stabilizing control law for the cascaded system that is comprised od the 
error dynamics (22) and the engine dynamics (9). To this end we assume that the engine dynamics is of full 
relative degree n, and therefore can be represented in the normal form without internal dynamic 

£ -1 ±e(0 = (A) - b 0 k T )x e (t) + b Q [i x>u e (t ) + d e (t)\ 

y e (t) = cqT x e (t) . (24) 

where the triplet (A, bo, Co) is the canonical form representation of chain of n integrators, k is an unknown 
constant vector, u> is the unknown constant control effectiveness of known sign (assume is positive), and 
y e (t) = A T(t) and u e (t) = A 

A. Prediction Models 

The prediction model for the lateral-directional tracking error dynamics is introduced according to equation 

e(t) = A m e(t) + B A(t)u(t) + Kj (t)x(t) + Kj (t)r(t) + W T (t)tp(t) + h(t) , (25) 

where the variables with ”hat” are the estimates of the corresponding unknown constant parameters, 
u(t) = [5 n (t) — 8®(t) AT(t) — AT°(<)] t , 8®(t) and A T°(t) are baseline aileron deflection and differen- 
tial thrust commands, A T(t) is the prediction of the differential thrust command to be design shortly, 
h(t) = [di(t)sign(xi(t)) CT 2 (i)sign(x2(i))] T , cr(t) = {d\ (t) is the estimate of the unknown bound er* 

and x{t) = [xi(0 X 2 (i)] T = B T Pe(t), and e(t) = e(t) — e(t) is the prediction error. The estimates A (t), 
K x (t), K r {t ), W T (t) and &(t) are given by the adaptive laws 

a (t) = 7in|A(t),M(t)x T (t)| 

K x {t) = 72 x(t)x T {t) 

K r (t) = 73 r(t)x T {t) 

W(t) = 74 V>(t)x T {t) 

&r(t) = lb\x{t)\, (26) 

where 71 > 0, 72 > 0, 73 > 0, 74 > 0 and 75 > 0 are the adaptation rates, the projection operator 17 
II (-, •) is introduced to keep the estimate A (t) positive, and the absolute value of the vector x(t) is defined 
componentwise. The prediction error dynamics can be written as 

e(t) = A m e(t) + B A(t)u(t) + A(t)u(t) + Kj(t)x(t) + Kj (t)r(t) + W T (t)ip(t) + a(t) - h(t) , (27) 

where A (t) — A — A (t), K x (t ) = K x — K x (t ), K r {t) = K r — K r (t), and W(t) = W — W(t) are the parameter 
estimation errors, u(t) = [0 A T(t)] T and A T(t) = A T(t) — A T(t) is the differential thrust prediction error. 

To be able to introduce a suitable prediction model for the engine dynamics in (24), we introduce filters 
following the conventional linear system parametrization scheme (see for example 15,18 ) 

e _1 €i(i) = (A -6°feo)^iW +b°y a (t), £i(0) = 0 

£ -1 £ 2 ( f ) = {A 0 -b 0 ko)€ 2 (t) + b°u e (t), ^ 2 (0) = 0 
e~% (t) = (A 0 -b°k^)^(t)+b°d e (t), € 3 (0)=0, 


( 28 ) 



where b° = [1 6° ... 6°] T is chosen such that the pair (Ao,b°) is controllable and the polynomial s n 1 + 

b 0 s n~2 + ... + b 0 n 

is Hurwitz, fco G K™ is chosen such that Aq — b°kQ is Hurwitz. Then, the differential thrust 
can be represented as 

Veit) = q T £ 1 (t) + + a e {t ) , ( 29 ) 

where we denote q = ko — k and a e (t) = £31 (t). Since the filters are stable, from the boundedness of the 
inputs the boundedness of the outputs follow. That is 07(f) is bounded by some positive constant a) since 
d e (t) is assumed to be bounded. The prediction model for the differential thrust is introduced as 

Ve{t) = q T (i)£i W + u>(t)€ 2 i{t) + d- e (t)sign(x 2 (i)) - ( 30 ) 

where q(t) and u>(t) are the estimates of unknown parameters q and u> respectively, and 07(f) is the estimate 
of the unknown bound a). The prediction error y e (t) = y e {t ) — y e (t ) is given by the equation 

Veit) = q T (f)£ i(f) + oj(t.)&i(t) + a e (t) - <r e (t)sign(x 2 (t)) , ( 31 ) 

where q(t) = q — q(t), and Q(t) = to — w(t) are the parameter estimation errors. The adaptive laws for the 
estimates q(t) and u>(f) in the differential thrust prediction model are defined as 

q(t) = 76 [j/e(*)£l(*) +A 2 (f)X2(i)£l(t)] 

£j{t) = 77II |w(f), y e (f)£ 2 i(f) + A 2 (t)x 2 (t)^2i(i)} 

d e {t) = M\Ve{t)\ + A 2 (t)|x 2 (i)|] , ( 32 ) 

where je > 0 , 77 > 0 and 78 > 0 are the adaptation rates, and the projection operator 17 II (•, •) is introduced 
to keep Cj(t) positive. 

The following lemma guarantees the closeness of the prediction models to the corresponding dynamics. 

Lemma 1 . The adaptive laws ( 26 ) and ( 32 ) guarantee boundedness of the estimates A (t), K x (t), K r (t), 
W T (t) and <x(f) as well as q(f), cD(f), 07(f) and the prediction errors e(f) and y e {t). Moreover, e(f) £ £ 2 
and y e (t ) £ £ 2 . 

Proof. Consider the following candidate Lyapunov function 

V(t) = e T (t)Pe r (t) +tr|7f 1 A T (f)A(f) + "/f 1 Kj (t)K x (t) + ^Kj (t)K r (t) + j^ 1 w T (t)W(t)^ 

+ 7 5 " 1 ( <t * - °r{t)) T {cr* r - dr r (f)) +7 6 - 1 qT(f)q(f) +7f 1 w(t)w(t) -07(f)) 2 . ( 33 ) 

The derivative of V(t) is computed along the trajectories of systems ( 27 ), and ( 26 ), ( 31 ) and ( 32 ). 

V(t) = —ej (; t)Q r e r (t ) + 2 x T (f) (A(f)'u(t) + A (t)u(t) + Kj (t)x(t) + Kj ( t)r(t ) + W T (t)tp(t) + a(t) 

- h(t)] +2tr|7 1 ” 1 A T (t)A(t) + {t)K x (t) + (t)K r (t) + jf) 1 w T (t)W{t)} 

- 27s- 1 (<r* - &(t)) T cr{t) - 27 ^q 1 (t)q{t) - 2 7f 1 w(f)u>(f) - 2^ 1 (a* e - a e (t))a e (t) 

= -e T (f)Qe(f) + 2tr{A T (f) w(f)x T (f) ~ 7i _1 A(t) |+2tr|^J(f) x(t)x T {t)-lf l K x (t) j 

+ 2tr r(t)x T {t)l3 l K r (t) }+2tr{lT T (f) <p(t.)x T (t) - (t) } 

+ 2x 2 (t)A 2 (t) q T (t)£i(t) + u)(t)^2i(t) + — <j e (t)sign(x 2 (t)) + 2;x(t) T (<x(f) - h(t)) 

- 27 5 _1 (£r* - &(t)) T cr{t) - 2 ^ 1 q T (t)q{t) - 2 jf 1 Cj{t)u(t) - 2 % 1 (a* e - a e (t))a e (t ) . ( 34 ) 

where we use the well known trace property tr(a T b) = tr(ba T ) for any vectors a and b of the same dimension. 
Substituting the adaptive laws and after some algebra we obtain 

V(t) < —e T (t)Qe(t) + 2x(f) T (er(f) - h(t)) - 2(er* - &(t.)) T \x(t)\ - 2 y e (t) [q T (t)£i(t) + w(t)£ 2 i(t)] 

+ 2x 2 (f)A 2 (f)[o- e (t) - <T e (t)sign(x 2 (i))] - 2 (a* - b e (t))[\y e (t)\ + A 2 (f)|x 2 (i)|] , ( 35 ) 



where we have used the following properties of the projection operator 1 ' 

tr{A T (t) u(t)x T (t) -n|A(t),M(t)x T (t)} } <0 

u(t) y e {t)&i(t) + \ 2 (t)x 2 (t)&i(t) -n|(h(t), y e (t)&i(t) + A 2 (t)x 2 (<) 6 iW} <0, (36) 

The terms x(f) T (er(f) — h(t)) — (cr* — <x(t)) T |x(f)| are evaluated componentwise as follows 

X(f) T (<x(f) - h(t)) - (cr* - d-(t)) T \x(t)\ = Xi(t)cri(t) - \xi(t)\vi(t) - K - cti ( i))|Xi(*)l + X 2 (t)cr 2 (t) 
— 1x2 (i)|o- 2 (t) - (o -2 -°2(t)) |x 2 (*)l =Xi(t)vi(t) -°l\Xi(t)\ +X2(t)a 2 (t) - a* 2 1 % 2 (t,)\x2 (t) <0. (37) 

Taking into account the positivity of the estimate A 2 (t) the last two terms are evaluated as follows 

X 2 (t)X 2 (t)[a e (t) - d e (t)sign(x 2 (t))] - (a* - a e (t))[\y e (t)\ + A 2 (t)|X 2 (*)|] = cr e (t)x 2 (t)A 2 (t) 
-c r elX 2 (i)|A 2 (t) - (cr* - cr e (t))\y e (i)\ < -(a* - d e (t))\y e (t)\ < -(cr* - a e (t))y e (t ) . (38) 

Taking into account the inequalities above and the error definition (31) we obtain 

V(t) < -e T (t)Qe(t) - f a (t) , (39) 

which implies that the error signals A(f), K x (t), K r (t), W(t), cr* — ef(f), q(t), u> (t) and a* — cr e (t) are 
bounded. Then, the estimates A(f), K x (t), K r (t), W(t), a (t), q(t), cu(t) and a e (t) are bounded as well. At 
this moment we cannot conclude anything about the boundedness of e(t), e(t), y e (t), y e (t ) or y e (t ). However, 
we can integrate the inequality (39) and obtain 

f e T (r)Qe(r) + y 2 e (r) dr < V(0) - V(t) . (40) 

Jo L 

From. Lemma 1 it follows that V(t ) is bounded, therefore e(t) G £ 2 and y e (t ) G £ 2 . □ 

We notice that Lennna 1 holds independent of specific control design, which we provide in the next 
subsection following the steps from Reference . 19 

B. Control Signal 

Our goal is design aileron deflection and differential throttle commands S a (t) and AdV(t) for the cascaded 
prediction system 

e(t) = A m e(t) + B A(t,)u(t) + kj (t)x(t) + I<J (t)r(t) + W T (t)<p(t) + h(t) 

£ _1 ^ 2 (t) = (A 0 -b 0 kg)^. 2 (t) + b 0 u e (t), (41) 

such that e(t) — > 0 as t — > 00 . We notice that the second system in (41) is of relative degree one and 
minimum phase, and the matrix A$ — b^k^ is Hurwitz. Therefore the vector fco can be chosen such that the 
transfer function from the input u e (t) to the output £21 (t) has the form 

61 (*) = -^-Ue(s). (42) 

s + eu 2 

where uq > 0 and w 2 7 ^ 0. The rest of the system is input to state stable, where the input is £21 (t), and 
is bounded as long as £21 (t) is bounded. Therefor the control design besides the tracking task must also 
guarantee the boundedness of £21 (t)- The cascaded system under the consideration takes the form 

e(t) = A m e(t) + B A(t)u(t) + Kj ( t)x(t ) + Kj ( t)r(t ) + W T (t)tp(t) + h(t) 

£ - 1 6 i (t) = -W 1621 (t)+U2U e (t), (43) 

Since A (t) is diagonal with positive entries, we can design the control signal u(t) as follows 

u(t)=A~ 1 (t) -k](t)x(t)-Kj(t)r(t)-W T (t)ip(t)-h(t.) . 


(44) 



Whereas the signal 


S a {t) = «i(f) + S°(t ) , 


(45) 


where S°(t) is the baseline signal from (16) is the actual aileron deflection command, the signal A T(t) = 
u 2 (t) + A T°(t), where A T°(t) is the baseline signal defined in (16), represents differential thrust command 
to the engine dynamics. Since we will design the differential throttle command from the point of view of 
controlling the model (52), the stabilizing function is designed as follows 


MO 


u 2 {t) - g T (t)Xi (0 - < 7 e (t)sign(x 2 (Q) 

MO 


(46) 


which is well defined since c D(t) > 0 is guaranteed by the adaptive law (32). Next we filter the stabilizing 
function fii (t) through a second order stable filter 

£ _1 ii(t) = UfZ 2 (t) 

Z 2 (t) = -2Quj f z 2 (t)-uj f [z 1 (t)- (47) 

and introduce an error variable 77 (f) = £21 (0 — Z\ (f ) . The tracking error prediction model now can be written 
as 


e(0 = A m e(t ) + b 2 \ 2 (t)u>(t)r](t) + b 2 X 2 (t)ui(t)[zi(t) - n(t)] , (48) 

where b 2 is the second column of matrix B. Filtering the auxiliary error term X 2 (t)Lo(t)[zi(t) — /r(f)] through 
the stable filter 


C (0 = A m C{t) + b 2 Cj{t)r}i(t)[zi{t) - /i(t)] , (49) 

we introduce the compensated error e(f) = e(f) — £(t) which satisfies the dynamic equation 

e(t) = A m e(t) + b 2 \ 2 (t)u{t)r]{t) . (50) 

Differentiating 77 (f) we obtain 

£~ lr q(t) = -£ -1 wi$21 it) + u 2 u e (t) - LOfZ 2 (t ) , (51) 

which together with the compensated error dynamics (49) make a singular perturbation problem 13 

e{t) = A m e(t) + b 2 X 2 (t)u>(t)ij(t) 

= -wi6i(i) + u 2 u e (t) - u f z 2 {t) . (52) 

In Reference 19 we prove that the singular two time scale separation principle can be applied to the cascaded 
system (52) and the control signal u e (t) can be designed through from the perspective of controlling the 
reduced system 

— = - u>fz 2 (t ) - u>i£ 2 i (t) + u> 2 u{t) . (53) 

CLT 

in the time scale r = et. The boundary layer system in this case has the form 


y(t) = A m y(t ) , 


(54) 


where y(t) = e(t) — A.^b 2 X 2 {i)Cd{t)y(1:) , and is exponentially stable. The control signal is designed as 


= (*) = 


-cy(t) + (ui - c)&i(t) +WfZ 2 (t) 


(55) 


and exponentially stabilizes reduced system (53). 


co 2 



Then it can be proved following the Reference 19 that the solutions of the cascaded system (52 and the 
following system 

e(t) = A m e(t ) + b 2 \ 2 (t)u>(t)r](t) + b 2 \ 2 (t)Cb(f)[zi(t) - p(t)\ 


£ 1 r t (t) = 

-Wl$21 (t) +i0 2 U e (t) - C0fZ 2 (t) . 

(56) 

satisfies the relationship 




nr 

0 

11 

-to 

(57) 

for all t > 0 , and the relationship 




e(t) = 0(e) + O(to f : ) 



nr 

0 

11 

nr 

(58) 


for all t > t* , where t* depends on the convergence rate of the compensated error system (52) in faster time 
scale, and on e. From the boundedness of e(t) and 77 (f) and from the inequality (40) it follows thate(f), 
y e {t) and y e (t ) are bounded. Then, from the dynamics (27) it follows that e(t) is bounded. Therefore, the 
application of Barbalat’s lemma results in e(t) — ► 0 as t — > 00 . From this convergence one can infer that the 
relationship 

e(t) = O^+OfuJ 1 ) (59) 

holds as t — > 00 . 

Since the lateral-directional tracking error □ (t) is bounded, from the boundedness of the reference com- 
mand it follows the the lateral-directional variables are bounded, and the bounds are independent of 

the unknown parameters. That is the compact set f l XD can be chosen independent of the adaptive scheme 
such that xd(T) £ f l XD for all t > 0. From this considerations and Remark 3 follows the validity of the 
approximation (21). This completes the proof of the following theorem. 

Theorem 1 . Let the outer loop variables be in the compact set Ll XL . Then the differential throttle 
command ASr(t) given by the equation (55), along with the aileron deflection command S a (f) (45), adaptive 
laws (26), (32), stabilizing function (46) and filters (47), (49), (50) regulates the error signal e^ft) to 
a neighborhood of the origin with the diameter proportional to e and inversely proportional to the filters’ 
frequency to / . 


VI. Simulation Results 


For the simulation we consider a damaged GTM model that before the damage was in straight and level 
flight at h = 30000 ft with speed M = 0.8. We consider only linearized lateral-directional dynamics of the 
vertical tail loss GTM model, since the longitudinal dynamics do not change much and the model ia still 
symmetric about the vertical plane. The linearized lateral directional system is given by the matrices 


' -0.0794 

-0.9975 

0.0496 

0.0404 " 


0.0024 

0.0000 ' 

-0.5646 

-0.0955 

-0.0217 

0 

> Bd = 

-0.0523 

0.0293 

-7.4175 

0.2784 

-1.7186 

0 

6.4318 

0.0010 

0 

0.0507 

1.0000 

0 


0 

0 


It is unstable with the poles at (—1.6473, — 0.7007, 0.3683, 0.0862). The engines are given by the linear 
system 


0 

1 

0 

0 


0 

0 

0 

1 

0 

5 b e — 

0 

0 

0 

0 

1 


0 

_ -150.0135 

-171.4720 

-73.4960 

-14.0000 _ 


_ 150.0135 _ 


(61) 


which correspond to the normalized thrust and normalized throttle input. The engine dynamic model has a 
settling time of t B = 10.3851sec. 



We run a simulation assuming that the systems are known, but engine states are not accessible and only 
thrust output T n {t) is available for feedback. The engine dynamics are written in the equivalent form 

ii(i) = z 2 (t) 

*2 (t) = z 3 (t) 

Z 3 (t) = Z 4 (t) 

z 4 {t) = -0.586021 (i) - 2.6793 z 2 (t) - 4.59352 3 (» - 3.5 z 4 {t) + 0.5860it„(f) 

T n (t) = Zl (t), (62) 

for which an observer is constructed with the poles at (—8.9, — 8.8, — 8.7, — 8.6), and the states of the 
observer are used in generating the differential thrust as well as the aileron deflection commands. 

In the first simulation we stabilize the airplane given from initial conditions /3(0) = 1.5°, r(0) = 0,p(o) = 
0, 0(0) = 0. The simulation results are presented in Figures 3(a), 3(b), 4(a) and 4(b). 






(a) Sideslip angle and yaw rate 


(b) Roll angle and roll rate 


Figure 3. Lateral directional performance of the airplane 



Engine Response 





(a) Differential Thrust and Aileron deflection commands 


(b) Engine model performance 


Figure 4. Stabilizing control signals 


It is evident from the figures that the stabilization is achieved with the chosen engine model. The initial 
rapid reduction in sideslip angle is due to aileron deflection at the initial period as can be seen from Figure 
4(a). This also results in the negative banking of the airplane (Figure 3(b)), which buys some time for the 
engines to generate a differential thrust. However, the thrust generated it is still not sufficient enough to 


decrease the sideslip angle. When the aileron returns to a small region around equilibrium, the sideslip angle 
increases again for the next 10 sec, but the engines are able to generate enough torque to stop the increase 
and eventually regulate the sideslip angle. 

The differential thrust command is large during the initial period and results in the saturation of throttle 
inputs for both engines as can be seen in Figure 4(b). The engine responses are slow, but sufficient enough 
to compensate for the chosen initial disturbance. It should be mentioned that the airplane departs from the 
course during the recovery. After a stabile flight regime is established, the course can be changed, which we 
demonstrate in the next simulation. 

The next simulation results show the tracking performance of the airplane with initial —0.5° sideslip 
angle following the 15° bank angle command while regulating the sideslip. The performance is displayed in 
Figures 5(a), 5(b), 6(a) and 6(b). 


Dutch Roll Mode Roll Mode 



(a) Sideslip angle and yaw rate 


(b) Roll angle and roll rate 


Figure 5. Tracking performance of the airplane 



(a) Differential Thrust and Aileron deflection commands 


(b) Engine model performance 


Figure 6. Tracking control signals 


It can be seen from the figures that the proposed control scheme has also tracking capabilities. The 
steady turn rate is about 0.6° per second, which is sufficient enough to change the course of the given flight 
conditions as required. The required differential thrust is again higher in the initial period and results in the 
throttle saturation of both engines 4(b). 

We run the last simulation to stabilize the airplane with the proposed adaptive control when the engine 
dynamics are assumed to be unknown, but the relative degree is known. The results for the initial conditions 
/ 3 ( 0) = 1.5°, r(0) = 0 ,p(o) = 0,^(0) = 0 are displayed in Figures 5(a), 5(b), 6(a) and 6(b). Although the 
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Dutch Roll Mode 


Roll Mode 






(a) Sideslip angle and yaw rate (b) Roll angle and roll rate 

Figure 7. Adaptive stabilization of the vertical tail loss airplane 



Engine Response 





(a) Differential Thrust and Aileron deflection commands 


(b) Engine model performance 


Figure 8. Adaptive control signals 


convergence is somewhat slower than in the known case, the overall performance is similar and the airplane 
behaves in a similar way. 


VII. Conclusion 

We considered a control problem for the vertical tail loss generic transport aircraft model, using the en- 
gines’ differential thrust in place of a rudder, which is assumed to be lost. Although the required differential 
trust generation is substantially slower than the conventional yaw moment generation by the rudder deflec- 
tion, we show that the aircraft can still be locally controlled, using the proposed control algorithm. This 
approach is demonstrated in the simulations for both stabilization of the damaged airplane and tracking of a 
given command. Future validation of the approach will be conducted using full nonlinear aircraft simulations 
with realistic and accurate engine model. 
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